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Using a combination of computational simulations, atomic-scale resolution imaging and phe¬ 
nomenological modelling, we examine the underlying mechanism for nanodomain restructuring in 
lead zirconate titanate (PZT) nanodisks driven by electron beams. The observed subhertz nan¬ 
odomain dynamics are identified with relaxation oscillations where the charging/discharging cycle 
time is determined by saturation of charge traps and nanodomain wall creep. These results are 
unusual in that they indicate very slow athermal dynamics in nanoscale systems. 


The importance of finite-size boundaries for static po¬ 
larization configurations has been emphasized in flux- 
closure [1-5] and faceted [6-8] domain pattern studies 
in several ferroelectric materials. Recently it has be¬ 
come possible to reorient spatially domains and defects 
on nanoscales with controlled high resolution transmis¬ 
sion electron microscopy (HRTEM) [9-11]. Nanodomain 
dynamics in ferroelectric nanodisks can also be probed 
with real-time HRTEM, and have been observed to be 
considerably slower (10 sec) [11] than those in their nano- 
magnetic counterparts (50 nanoseconds) [12]. Although 
in principle, such slow nanodomain reorientation could 
arise from thermal effects, it has been shown that such 
heating is negligible for our experimental conditions [13]. 
Here we emphasize that the observed driven domain dy¬ 
namics of such free-standing ferroelectric nanostructures 
are strongly influenced by electrostatic conditions. In a 
combined computational, experimental and phenomeno¬ 
logical effort, we identify the charging/discharging mech¬ 
anism underlying the observed oscillations between circu¬ 
lar and hexagonal nanodomain patterns in lead zirconate 
titanate (PZT) nanodisks driven by electron beams. Us¬ 
ing phase-field modelling [10, 14, 15], we demonstrate the 
crucial influence of charge on the nanodomain patterns, 
and show that we can reproduce domain structures sim¬ 
ilar to those observed in experiment. The subhertz fre¬ 
quency response of the nanodomain reorientation results 
from relaxation oscillations [16-18] between threshhold 
charging of surface traps [19] and subsequent discharg¬ 
ing by nanodomain wall creep [20], where quantitative 


agreement is made using prior measurements on PZT 
samples. Experimentally edge-supported 8-nm diameter 
disks of PZT were exposed to constant HRTEM beams 
with a 0.5 A/m^ probe current density where further 
detail can be found elsewhere [11]. Two distinct nan¬ 
odomain patterns were observed: the nanodomain walls 
“flop” between being normal to the disk perimeter to be¬ 
ing parallel to it as displayed in Figure la and Ic. Since 
HRTEM only reveals ferroelastic domain walls, modelling 
is needed to examine the full charging role played by the 
electron beam that drives the observed faceting oscilla¬ 
tions and the simultaneous realignment of nanodomain 
walls. It is known that electron beams can lead to charg¬ 
ing of insulating materials [21]. Eurthermore in the con¬ 
text of ferroelectrics, it has been shown that electron 
beams can be used to switch the polarization direction 
[22] and can even stabilize unusual quadrant patterns in 
free standing nanodisks [10]. In order to investigate e- 
beam induced charging effects in ferroelectric nanodisks, 
we use a phase-field approach to simulate how domain 
patterns are influenced by radial fields generated by uni¬ 
form free charge. The model, based on earlier studies 
[10, 15, 23], has been adapted to incorporate the circular 
geometry (see supplementary material [24]) of the ex¬ 
periment; the goal is to determine whether the observed 
nanodomain pattern realignment is driven by depolar¬ 
ization fields in a charging/discharging cycle. We first 
simulate the uncharged case by solving the phase field 
equations for a system quenched from the paraelectric 
state. A uniform charge density, A^e = 4 x 10^® m“^, is 
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FIG. 1. HRTEM images for (a) circular (unfaceted) and (c) 
faceted 8-nm PZT dots; the ferroelastic domain walls are il¬ 
lustrated as purple lines, and the angles between them are 
indicated. Results from simulations (distribution of polariza¬ 
tion vectors) for (b) uncharged and (d) charged nanodisks; red 
lines indicate the approximate positions of domain walls. Do¬ 
main walls are perpendicular to the surface in the unfaceted 
cases (a. and b.), but parallel to the surface in the faceted 
cases (c. and d.) 

then applied and the evolution of the domains is moni¬ 
tored. Fig 1(b) shows the simulated pattern for the un¬ 
charged case with a flux closure domain pattern where 
the solid red lines indicate the positions of the ferroelas¬ 
tic domain walls, normal to the disk perimeter. We note 
the excellent agreement with the domain wall configura¬ 
tion in the HRTEM image in Fig. 1(a), suggesting that 
it corresponds to an uncharged disk. Figure 1(d) shows 
the nanodomain pattern when the nanocrystal is charged 
{Ne = 4 X 10^® ni“^), i.e. during the charging process, 
the flux closure domain pattern transforms to one where 
polarization vectors have a non-zero radial component at 
the perimeter and a polarization vortex type state, char¬ 
acterized by a vanishing radial component, in the center. 
This may be understood in terms of the radial electric 
field that vanishes at the center. This radial field arises 
due to the uniformly charge density [10] associated with 
the electron beam induced charging. Vortex formation 
at the center occurs to avoid head-to-head polarization, 
resulting in areas having P • n = 0 (center) and P • n ^ 0 
along the perimeter. Here, a radial unit vector is intro¬ 
duced as n = |[[~[]°| , where Tq = {xo,yo) represents the 
coordinates of the centre of the disc. Regions with such 
distinct P • n are separated by ferroelastic nanodomain 
walls that are accessible to the electron beam technique. 
The solid red lines in Fig 1(d) indicate such ferroeleastic 
domain walls which are aligned along the perimeter of 


the nanocrystal. Thus this situation corresponds to the 
HRTEM image in 1(c), indicating that it results from 
nanodisks with charged boundaries. 

While this charging/discharging scenario describes the 
realignment of domain walls, it does not directly model 
the observed kinetic faceting. We emphasize its dynam¬ 
ical nature since two-dimensional faceting cannot occur 
in thermal equilibrium at finite temperatures [25]. How¬ 
ever, in Fig 1(d), we see that the distribution of P • n 
is anisotropic, namely strongly 0-dependent, and this is 
further emphasized graphically in Figure 2. 
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FIG. 2. Distribution of the normalized radial component of 
the polarization, ;^P • n corresponding to Fig Id. The ra¬ 
dial component of the polarization is zero in the center and 
nonzero towards the edges. Note that the distribution shows 
strong anisotropy (i.e. 0-dependence) 

Indeed the distribution of elastic energy densities will also 
be orientation-dependent due to the underlying crystal 
axes of PZT. More generally kinetic faceting is dynami¬ 
cally selected to accommodate external anisotropies that 
could also include applied stress and substrate-related 
issues [6-8] where details depend on the specifics of the 
experimental situation and the materials involved; in¬ 
deed faceting nanodomain patterns were not observed in 
HRTEM experiments on PZT nanotargets with square 
and triangular boundaries [11]. We emphasize that nan¬ 
odomain realignment and faceting is observable only be¬ 
cause our experiments uniquely involve free-standing, 
edge-supported nanodisks [11], since substrate interac¬ 
tions would prevent such domain restructuring. Our 
phase-field simulations indicate that this nanodoman re¬ 
alignment occurs between charged and uncharged states 
of the PZT nanodisk. When the disk is uncharged, the 
usual flux closure is observed implying that the ferroe¬ 
lastic domain walls are normal to the surface. When it 
is charged, the polarization at the surface becomes ra- 
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dial and a vortex is formed at the center. In HRTEM, 
such a state is revealed by its nanodomains parallel to the 
surface. The faceting is selected to minimize both polar¬ 
ization and external anisotropies, principly due to strain 
fields at the perimeter due to the underlying (111) crystal 
axes of PZT. This electrostatic analysis accounts for the 
two observed nanodomain configurations, but suggests 
a frequency response on gigahertz scales in contrast to 
the subhertz frequencies observed [11]. Because thermal 
effects have been ruled out for these experimental con¬ 
ditions [13], this slow time-scale suggests extreme non¬ 
linearity and possibly mechanical overshoot. Our phase- 
field modelling indicates a driven and reversible charg¬ 
ing/discharging cycle that is naturally described as an 
effective relaxation oscillation circuit [17, 18]. This is a 
system that is intrinsically out-of-equilibrium: it cycles 
between an unstable state and a frustrated attempt to 
reach equilibrium that is prevented by onset behavior of 
its nonlinear element that forces it back to its original 
unstable state. Electronically the nonlinear circuit el¬ 
ement is often a neon light or a unijunction transitor, 
but its key feature is a threshhold voltage for negative 
resistance in parallel with a conducting channel that is 
usually a capacitor. 

To our knowledge the first report of self-sustaining sub¬ 
hertz oscillations in a physical system was in a wide-gap 
semiconductor ZnO [16], and we propose that the e-beam 
experiments on PZT nanodisks also present such an ex¬ 
ample. It is known that electron beam irradiation creates 
charge traps and defects in semiconducting targets. The 
time-scale for the charging cycle will then be set by the 
charge saturation of these surface states. This charge 
storage threshhold in PZT targets exposed to electron 
beams has been measured quantitatively to be 5C/m^. 
For the experiments of interest [11] with current density 
of 0.5A/m^, this then corresponds to a charging time 
of 


"^charging — 


5C/m^ 

0.5C/s-m^ 


= 10 sec 


( 1 ) 


that is the same order of magnitude as the observed 
faceting oscillations [11]. We note that if we assume 
that these traps are singly charged, then there are 
1 6x^10'-^° C/e = 3 X 10^®traps/m^ that is consistent with 
concentrations reported for oxygen vacancies in commer¬ 
cial quality spin-on PZT films [26]; this value is also in 
good agreement with the electron concentration used in 
our simulations {Ng. = 4 x 10^®m“^) for 8 nm disks of 
thickness 100 nm. Finally we note that our system is a 
charge analogue of a seesaw with a faucet continuously 
dripping water into a container on one end, and a block 
on the other one. In this textbook mechanical relax- 
ational oscillator [27], the seesaw flips when the weight 
of the water container is greater than that of the block; 
however once the container touches the ground it emp¬ 


ties, the seesaw goes back to its original position and the 
whole process starts again. In the context of our present 
work, the electron beam is analogous to the water flow 
from the faucet; the threshold is related to the Zener 
breakdown [28] that occurs in pn-junctions at low helds 
ca. E ^ l;j[^ and involves excitation of electrons to the 
conduction band from the valence band or trap states 
within the bandgap. It is non-destructive and reversible, 
and thus differs from avalanche breakdown that occurs 
at higher helds. The charging of the nanocrystal creates 
depolarization helds that realign the domains and lead 
to the observed faceting (Fig 1). During the initial part 
of the cycle, the current from the electron beam is hlling 
the surface states with charge and the nanodisk behaves 
like a capacitor. However once the charge saturation 
threshhold for these traps and defects in the PZT disk 
is reached [19], the depolarization helds become strong 
enough to initiate the current how in the nanodisk. This 
current leads to discharging and the subsequent decrease 
of the depolarization held that caused the faceting; con¬ 
sequently the domains realign back to their uncharged 
conhgurations. The presence of surface states and large 
electromechanical coupling means that this involves do¬ 
main creep in a random elastic medium that has been 
discussed quantitatively and measured in PZT hlms [20]. 
For the experimental parameters of interest, the domain 
wall creep in PZT epitaxial hlms is found to w = 10“® 
m/s [20] so that for the dots of roughly 10 nm (of thick¬ 
ness 100 nm), the discharging time-scale is 


^discharging — 


lOnm 

10“®m/s 


= 10 sec 


( 2 ) 


that is again in the same order of magnitude as the ob¬ 
served time-scale. This charging/discharging cycle is ex¬ 
pected to continue as long as the electron beam is on and 
a there is a continuous how of charge into the nanodisk. 
In contrast with many relaxation oscillators, in these 
PZT nanodisks there are slow time-scales associated with 
both the charging and the discharging cycles where the 
rate-limiting one is clearly the slower of the two. Thus 
the cycling times, achieved both by macroscopic (charg¬ 
ing threshhold) and microscopic (discharging via domain 
creep) methods, is consistent with measurement. We 
note here that the two underlying phenomena, charge 
trapping and slow domain creep due to electromechanical 
coupling and mass, are both absent in magnetism which 
explains the many orders of magnitudes that separate the 
frequency response of field-driven ferroelectric and ferro¬ 
magnetic nanodisks. In summary we have studied the 
mechanism underlying nanodomain restructuring in PZT 
nanodisks driven by electron beams. Our phase-field 
modelling identifies the two observed nanodomain pat¬ 
terns with uncharged and charged free-standing bound¬ 
aries; furthermore anisotropy in the charged case selects 
faceting at the perimeter to minimize electrostatic and 
elastic energy costs. The observed subhertz frequency 
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response, many orders of magnitude lower than expected 
from an electronic mechanism, is explained as relaxation 
oscillations where the charging/discharging times are de¬ 
termined by charge trap saturation and nanodomain wall 
creep. We predict that the oscillation time-scale can 
then be tuned by changing the surface trap density, ei¬ 
ther by varying the beam current or by annealing the 
PZT nanodisks in oxygen/ozone environments. Though 
this charging/discharging cycle here is driven by electron 
beams, similar behavior could be achieved by gating the 
PZT nanodisks to achieve a voltage-controlled source- 
gate-drain device. It would be unusual to have such small 
semiconductor triodes with such slow time-scales, and 
they could be important for biological applications like 
medical implants where time-scales are naturally of the 
order of seconds. 
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